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Abstract 

Polarized deep inelastic scattering (DIS) data are analyzed in leading and next-to- 
leading order of QCD within the common 'standard' scenario of polarized parton 
distributions with a flavor-symmetric light sea (antiquark) distribution Sq, and a 
completely SU(3) f broken 'valence' scenario with totally flavor-asymmetric light 
sea densities (5u 7^ 5d 7^ 5s). The latter flavor-broken light sea distributions are 
modelled with the help of a Pauli-blocking ansatz at the low radiative/dynamical 
input scales of ^lo(nlo) = (0-40) GeV 2 which complies with predictions of the 
chiral quark-soliton model and expectations based on the statistical parton model 
as well as with the corresponding, well established, flavor-broken unpolarized sea 
(d > u). Present semi-inclusive DIS data cannot yet uniquely discriminate between 
those two flavor-symmetric and flavor-broken polarized light sea scenarios. 



1 Introduction 

The polarized parton distributions of the nucleon have been intensively studied in recent 
years [1 - 14]. The conclusion has been that the experimental data dictate a negatively 
polarized antiquark component, and show a tendency toward a positive polarization of 
gluons. Presently we possess a lot of precise data [15 - 24] on the polarized structure 
functions of the nucleon, some of them very recent |23|, which justify a renewed 
investigation of the aforementioned issue. This alone, however, does not provide the main 
motivation for this project, rather the improved understanding in recent years of the 
situation in the unpolarized parton sector [25 - 28] provides important insights for the 
corresponding polarized parton densities. In particular, one notes that the unpolarized 
sea (antiquark) distributions are flavor-asymmetric (d > u), which can be understood in 
terms of flavor mass asymmetries and Pauli-blocking effects p9[ [3(J. The main objective 
of the present paper is to transcribe these insights into the polarized parton sector as will 
be described in Section 3. In Section 2 we shall, for completeness, present an analysis 
within the framework of the simplified SU(3)f symmetric 'standard' scenario in which the 
flavor-asymmetries in the polarized antiquark sector are neglected. This is done in view 
of the fact that in many situations these flavor asymmetries are unobservable as is the 
case for (most) presently available data which cannot provide any reliable information 
concerning this issue. 

Measurements of polarized deep-inelastic lepton nucleon scattering yield direct infor- 
mation [15 - 24] on the spin asymmetry 



N = p, n and d = (p + n)/2 where in the latter case we have used gf = (g{ + gi)[l — 
3u D /2}/2 with uo D = 0.058, R = F L /2xF 1 = (F 2 - 2xF 1 )/2xF 1 and subdominant contri- 
butions have, as usual, been neglected. 




(1.1) 
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We emphasize that, as in our original analysis [l], we compute both g\ and F\ entirely 
in leading-twist QCD. In particular, in order to obtain F\, we use the parton densities of 



GRV98 [25| along with LO (note that R = at leading order) or NLO coefficient functions 



for F 2 and R in (1.1). An alternative, frequently adopted [[], f|, ||, ||, [L0], [L3| approach is 
to take F 2 (x,Q 2 ) and R(x,Q 2 ) from experimental measurements, which is motivated by 
the fact that leading-twist calculations of F 2 (x,Q 2 ) and R(x,Q 2 ) do not agree very well 
with experimental determinations in the region of low Q 2 and W 2 = Q 2 (l — x)/x. These 
regions are affected by power-suppressed contributions and are therefore excluded from 
all unpolarized DIS analyses. The presently available data in the polarized case, however, 
do not allow to impose similarly 'safe' cuts (Q 2 > 4 GeV 2 , W 2 > 10 GeV 2 ) without losing 
too much information. On the other hand, the Q 2 range accessed so far in polarized DIS 
does not allow for extracting the magnitude and shape of power-suppressed contributions 
reliably. To study the issue further, we performed fits to the Ax data in both possible 
ways, i.e. with leading-twist calculations of F 2 and R as well as with their experimental 
results, admitting at the same time an 'effective higher-twist' contribution to A\ in terms 
of a factor (1 + A(x)/Q 2 ) with A(x) to be determined by the data. The outcome of this 
analysis was that A(x) is consistent with zero if we use leading-twist QCD for F 2 and R, 
but that it is sizeable and important in the fit if F 2 and R are taken from experiment. We 
take this as an indication that our preferred approach is more consistent and less liable to 
modifications by higher-twist terms. This view is also corroborated by the fact that the 
DIS A\ data show only a very mild Q 2 -dependence, even toward low values of Q 2 . The 
consistency of the polarized parton densities as extracted from DIS data at comparatively 
low values of Q 2 with measurements of other hard processes at higher scales can be studied 
soon at RHIC and perhaps in the future at a polarized ep collider. 

In NLO, gf(x,Q 2 ) is related to the polarized (anti)quark and gluon distributions 
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5f(x, Q ) = /+ — /_ in the following way: 



Q 2 ) = \ £ <5 2 ) + <5 2 ) + * i5qN + 5f) + l 5Ca * 5a]} 

(1.2) 

with the convolutions 5Cf*f being defined in the usual way. The MS coefficient functions 
SCf(x) can be also found in |l|], where all necessary ingredients for the Q 2 -evolution 
have been formulated as well. A similar expression holds for the unpolarized structure 
function F^(x,Q 2 ) with its spin-averaged parton distributions f(x,Q 2 ) = f + + /_ and 
the unpolarized Wilson coefficients can be found, for example, in [|3l|. The LO and 
NLO(MS) input scales, running coupling constants and parton distributions, employed in 
the positivity constraints \5f \ < f, will be adopted from GRV98 [ 25] . Furthermore we 



shall, as always, use the notation 5q p = 5q and q p = q, and neglect the marginal charm 
contribution to stemming from the subprocess j*g — > cc ||32|| . The charm contribution 
to Ff 1 is also small in the kinematic range covered by present polarization experiments. 

The total helicity of a specific parton f = u, u, d, d, s, s, g is given by the first (n = 1) 
moment 

Af(Q 2 ) = [ dx5f(x,Q 2 ). (1.3) 
Jo 

Thus, according to (1.2), 

Ti(Q 2 ) = j\x gi (x,Q 2 ) = \J2< [ A ?W 2 ) + A?(Q 2 )] (l - (1-4) 

since AC q = — 3CV/2 = —2 and AC g = 0. Therefore we have in general 

' x _ a s (Q 



r?- n (Q 

with the flavor-nonsinglet components 



+ S A * + 5 AE( ° !) 



7T 



;i.5) 



Aq 3 = Au + Au-Ad-Ad (1.6) 
Aq 8 = Au + Au + Ad + Ad-2(As + As) (1.7) 
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being conserved, i.e. (^-independent, and the flavor-singlet component is given by 



A£(Q 2 ) = i A( l(Q 2 ) + A <?(Q 2 )] = A?s + 3 [As(Q 2 ) + As(Q 2 )] . (1.8) 



These quantities will be subject to various constraints (derived from hyperon /3-decays) 
depending on the specific model scenarios under consideration to which we shall turn in 
the next two Sections. 

Finally, the fundamental helicity sum rule reads 



where L q+g refers to the total orbital contribution of all (anti) quarks and gluons to the 
spin of the proton. 

2 SU(3)y symmetric 'standard' (unbroken— sea) sce- 
nario 

As stated in the Introduction, present data do not provide sufficient information concern- 
ing the flavor-asymmetries of the polarized sea distributions. Thus present day studies 
must, as a first approximation, neglect this issue unless one is willing to adopt some models 
for the flavor-asymmetries as will be done in Section 3. Here we follow the procedure 
presented in GRSV95 [|TJ. The searched for polarized NLO (as well as LO) parton distri- 
butions 5f(x,Q 2 ), compatible with present data [15 - 24] on Ai(x,Q 2 ), are constrained 
by the 'standard' sum rules 



± = ±AZ(Q 2 ) + Ag(Q 2 ) + L q+g (Q 2 ) 




Aq 3 



F + D = g A = 1.2670 (35) 



(2.1) 



Aq 8 



3F-D = 0.58 ±0.15 



(2.2) 
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where the updated values for F and D have been taken into account and the error 



estimate in Eq. (2.2) is due to |34|. Thus, Eq. (1.5) becomes 

V2 



T p { n (Q 2 ) 



±^(F + D) + ^(3F - D) + l - (As(Q 2 ) + Ar s (Q 2 )) 



1 a s(Q 2 

71 J 

(2.3) 

i.e. one needs here a finite sizeable strange sea polarization As(Q 2 ) < in order to achieve 



the experimentally required reduction of the Ellis-Jaffe LO expectation [35 



r?, E j = ^(F + D) + ^(3F - D) ~ 0.186 . (2.4) 

Furthermore, in the 'standard' scenario one assumes an unbroken SU(3)/ symmetric sea, 

5q(x, Q 2 ) = 5u = 5u sea = 5d = Sd sea , = Ss = 5s . (2.5) 

For the determination of the NLO (LO) polarized parton distributions Sf(x, Q 2 ) we follow 
our original analysis [Q by relating the polarized input densities to the unpolarized ones, 
using some intuitive theoretical arguments [^] as guidelines. We employ the following 
ansatz for the LO and NLO (MS) polarized parton distributions at an input scale Q 2 = /i 2 



0: 



5u(x,n 2 ) = N u x a ^l-xf u u(x,fi 2 ) GRy 

5d(x^ 2 ) = N d x^(l-x)^d(x,fi 2 ) GRY 

Sq(x,fx 2 ) = NqX a «(l-xf«q(x,fj. 2 ) GRY 

6g(x^ 2 ) = N g x a *(l-x)^g(x^ 2 ) GRy (2.6) 



with the LO and NLO unpolarized input densities referring to the ones of GRV98 [IE] 
and q = (u + d)/2 should be considered as the reference light sea distribution for the 
'standard' unbroken-sea scenario in (2.5). The parameters of our optimal LO densities at 
/i 2 Q = 0.26 GeV 2 and the ones of the NLO(MS) densities at /^nlq = 0.40 GeV 2 are given 
in Table I. These optimal LO and NLO (MS) fits correspond to a x 2 per degree of freedom 
(Xdf) o f Xdf, lo — 0-84 and to Xdf.nlo — 0.81. The polarized gluon density in (2.6) 
is, as usual, rather weakly constrained by present data. Although a fully saturated (via 



the positivity constraint) gluon input 5g(x,fi 2 ) = ±g(x, /i 2 ) is disfavored, a less saturated 
5g(x,fi 2 ) = ±xg(x,fi 2 ) input or even a vanishing (purely dynamical) input 6g(x,fi 2 ) = 
are fully compatible with present data. The latter choice, however, seems to be unlikely 
in view of Sq(x, /i 2 ) 7^ 0. 

In Fig. 1 our NLO results are compared with the data on A±(x, Q 2 ) as well as with our 
old original NLO (MS) fit ffl. The differences between these two results are small, except 
perhaps for A™ in the large-x region. Our new LO fit is similar to the NLO one shown in 
Fig. 1 by the solid curves. The Q 2 -dependence of our LO and NLO 'standard' scenario 
fits at various fixed values of x is shown in Fig. 2 and compared with all recent data 
on gi(x,Q 2 ), including the most recent E155 proton data ||24j| . The main reason for our 



LO results being larger than the NLO ones in the small-x region is due to the vanishing 
of R p (x,Q 2 ) in LO in (1.1). The corresponding x-dependence of g^(x,Q 2 = 5 GeV 2 ) 
is shown in NLO in Fig. 3 where the expected extrapolations into the yet unmeasured 
small-x region down to x — 10~ 3 are shown as well. The solid curves refer to our 
optimal NLO fit (with the input given in Table I) and allowing our optimal total x 2 m 
Table I to vary by one unit, 5\ 2 = ±1, gives rise to the shaded areas due to different 
choices of the polarized gluon input at Q 2 = /iN LO in (2.6) such as 5g = ±xg, etc. 
In particular a vanishing polarized gluon input Sg(x, /^nlo) = 0' ^ s f° r the ti me being 
entirely compatible with all present data as shown by the dashed curves. On the other 
hand, fully saturated (via the positivity constraint) gluon inputs Sg(x, /Un LO ) = ±5 appear 
to be disfavored by present data as shown by the dotted curves in Fig. 3. It should be 
furthermore noted that the shaded bands in Fig. 3 contain polarized gluon densities which 
correspond to first moments Ag(Q 2 = 5 GeV 2 ) between —0.81 and 1.73, according to input 
moments between Ag(filf LO ) = —0.45 and 0.7, respectively, i.e. even negative total gluon 
polarizations are compatible with present data. The same results hold of course also 
in LO. Future dedicated polarized small-x measurements and upcoming determinations 
of Sg(x, Q 2 ) should be useful in removing such extrapolation ambiguities caused by our 
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present poor knowledge of 5g(x,Q 2 ). 

Our corresponding LO and NLO parton distributions at the respective input scales 
Q 2 = /i 2 nlo in Eq. (2.6) with the 'standard' scenario fit parameters given in Table I 
are shown in Fig. 4. The main differences between our new input densities and our 
old GRSV95 ones Jl| are somewhat harder 5d (due to the new neutron data) and Sg 
distributions although, as discussed above, the polarized gluon distribution in Fig. 4 is 
only slightly preferred by our 'optimal' fit to presently available data. The polarized input 
densities in Figs. 4(a) and 4(b) are compared with our reference unpolarized valence-like 
LO and NLO dynamical input densities of pi| which satisfy of course the positivity 
requirement \Sf \ < f as is obvious from Eq. (2.6). It should be nevertheless emphasized 
that the parameters, resulting from our rather general LO and NLO fits, are always such 
that these positivity conditions are automatically satisfied, i.e. there is practically no need 
to impose them separately. The distributions at Q 2 = 5 GeV 2 , as obtained from these 
LO and NLO inputs at Q 2 = fi 2 , are shown in Fig. 5 where they are also compared with 
our old NLO GRSV95 []J results. It should be noted that the substantially harder input 
polarized gluon density xSg(x, fi 2 ), in particular in LO, in Fig. 4(a) as compared to our 
old GRSV95 fit, causes the sea density x5q(x, Q 2 ) to oscillate (slightly) in the large- X 
region at Q 2 > fi 2 as shown in Fig. 5 |38| . 

Next let us turn to the first moments (total polarizations) Af(Q 2 ) of our polarized 
parton distributions, as defined in (1.3), and the resulting r^' n (Q 2 ) in (1.5). It should be 
recalled that, in contrast to the LO, the first moments of the NLO (anti)quark densities 
do renormalize in the MS scheme, i.e. are Q 2 dependent (see, e.g., Ref. Jl|), whereas the 
gluon polarization Ag(Q 2 ) renormalizes in both cases. Our LO 'standard' fit implies 

Am = 0.871, Ad = -0.396, Aq = As = As = -0.054 
Ag(nl Q ) = 0.190, A#(5GeV 2 ) = 0.684, A#(10GeV 2 ) = 0.802 (2.7) 
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which result in AS = 0.259 and 



r? = 0.151, r? = -0.061. (2.8) 

Our NLO results are summarized in Table II at some typical values of Q 2 . Both our 
LO and NLO results for Y p { n (Q 2 ) are in satisfactory agreement with recent experimental 
determinations [17 - 24]. Furthermore, due to the constraint (2.1), the Bjorken sum rule 
|39| holds manifestly in LO and, according to (1.5), the NLO a s -corrected sum rule reads 

r?(Q 2 ) - r? (Q 2 ) = \g A {i- ^p-) ■ (2.9) 

It is also interesting to observe that at our low input scales Q 2 = nlo = 0.26, 0.40 
GeV 2 the nucleon's spin carried by the total helicities of quarks and gluons amounts only 
to 

±AZ + Ag^ 2 LO ) ~ 0.32 
l - AE ( /4 LO ) + Ag(^ h0 ) ~ 0.35 (2.10) 

which implies for the helicity sum rule (1.9) already a sizeable orbital contribution 
L q+g (fi1 Q NLO ) ~ 0.18, 0.15 at the low input scales. Although this is in contrast to our 
somewhat more intuitive previous GRSV95 result 0, L g+9 (/iLo nlo) — 0' ^ should be 
kept in mind that, for the time being, Ag(Q 2 ) is rather weakly constrained by present 
data as was discussed above. 

Finally, for completeness we have also performed a NLO analysis in a different factor- 



ization scheme, the so-called chirally invariant (CI) or JET scheme ||, [40], [44J], but any 
other choice would do as well for studying the scheme dependence of our MS fit results. 
Here, among other things P, [IT], O], the total helicity of quarks, AE CI , is conserved, i.e. 
Q 2 -independent, and is related to our AE in the MS scheme via 

AE(Q 2 ) = AE CI - 3 MQ 2 )ci • (2.11) 

Z7T 
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Similarly agreeable fits as the ones in Figs. 1-3 can be obtained in this scheme, e.g. by 
choosing a large positive gluon density with a total (input) helicity Ag(fj,^ LO )ci — 0.6 — 0.7 
being about three times larger than in Table II and the sea density 5qci(Aqci) turns out 
to be roughly 50% smaller than the one of our best fit in the MS 'standard' scenario; here 
the total quark helicity increases to ASci — 0.4. 

3 SU(3)/ broken 'valence' (broken— sea) scenario 

The assumption of the flavor symmetric 'standard' scenario with its unbroken sea density 
in (2.5) is expected to be unrealistic, following our experience in the unpolarized case 
where a suppression of the strange sea component is required, as accomplished by the 
vanishing input s(x,fi 2 ) = s(x,fi 2 ) = in GRV98 f2~5| , in order to comply with experi- 



mental indications f|2], [E| of an SU(3) / broken sea, and the positivity constraint 5s < s. 
Thus, in GRSV95 [0] we also considered a 'valence' scenario where, in contrast to (2.5), 

Sq(x, /i 2 ) = 5u = 5d = Su sea = 5d sea 

5s(x, fj 2 ) = Ss(x,ix 2 ) = 0. (3.1) 

Furthermore the full SU(3)j flavor symmetry, giving rise to the constraints (2.1) and (2.2), 
is broken in the 'valence' scenario to the extent |43J that the flavor-changing hyperon (3- 
decay data fix only the total helicity of valence quarks Aq v = Aq — Aq: 

Au v (fi 2 ) - Advi/j 2 ) = F + D (3.2) 
Au v {fx 2 ) + Ad v (fx 2 ) = 3F-D, (3.3) 

i.e. Aq 3 = Au v — Ad v and Aq s = 3F—D+4Aq at Q 2 = fi 2 [[!]] according to (3.1). Therefore 
a light polarized sea Aq < suffices here to account for the reduction of the Ellis-Jaffe 
estimate (2.4). This is the reason for our simplifying assumption of a maximally broken 
SU(3)/ strange sea input in (3.1) in order to reduce the number of input distributions 
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to be fitted to the rather scarce available polarization data, which are now sufficient for 
fixing these input distributions. (Future high-statistics data should allow, at least in 
principle, to extract the total strange sea polarization without employing any simplifying 
assumption, as for example from (1.8), A(s + s) = (AS — Ag 8 )/3 .) The quality of the 
fits obtained is comparable to that of the 'standard' flavor-symmetric scenario discussed 
and presented in the previous Section, cf. Fig. 1. We refrain, however, from presenting 
these results here explicitly because the assumed remaining flavor-isospin symmetry of 
the light sea components in (3.1) appears to be somewhat artificial and unnatural in view 
of the flavor-asymmetric unpolarized light sea distributions d(x, Q 2 ) > u(x, Q 2 ) [25 - 28]. 

Turning now to the presumably more realistic scenario where also the flavor-isospin 
symmetry of the polarized sea is broken, we note, as already pointed out in the Introduc- 
tion, that some model assumptions are needed for the corresponding input distributions. 
The analysis of the unpolarized structure functions yields 

d(x, /i 2 )/-u(x, jj 2 ) ~ u(x, fi 2 )/d(x, /i 2 ), (3.4) 



which holds to a rather good accuracy for the GRV98 distributions [p5 |. This propor- 
tionality relation is expected to hold approximately at least for 0.01 < x < 0.3 where 
the breaking of the light sea d > u is directly tested experimentally via Drell-Yan dilep- 



ton production in pp and pd collisions ||5|] and semi-inclusive 7r ± production in ep and 



ed reactions [[|6|]. The relation (3.4) may be considered [3C] as a manifestation of the 



Pauli-blocking effect [47| which should be relevant also in the polarized parton sector. 
We therefore estimate the flavor-symmetry breaking of the polarized sea to be given in a 
first approximation by [p0|1 

5d(x, jj 2 ) /5u(x, ji 2 ) = 5u(x, /j, 2 )/5d(x, fi 2 ) (3.5) 



together with the previously advocated 

Ss(x, jj 2 ) = 5s(x,n 2 ) = 
10 



(3.6) 



According to the unpolarized case above, we expect the proportionality relation (3.5) to 
hold approximately at least for 0.01 < x < 0.3. 

It should be reemphasized that, in complete analogy to unpolarized DIS, data on 
inclusive polarized DIS in kinematical regimes where only photon exchange is relevant 
give information only on the sums of quark and antiquark polarizations for each flavor, 
Sq(x, Q 2 ) + Sq(x, Q 2 ), appearing in g\(x, Q 2 ) in (1.2). This implies, in particular, that the 
amount of flavor-SU(2) breaking in the polarized sea cannot be determined from such data 
since one can always change the parton densities by Sq — ► 5q — 5(fi q and 5q — > 5q + 5<f) q for 
any arbitrary functions 8(p u (x) and 5<pd(x) without affecting at all the measured structure 
functions gi' n,d (x , Q 2 ) . Restrictions on 5(fi q (x) occur only in our 'valence' scenario where 
A(f>d = — A<j) u for the first moments, as in general implied by our modified (Au ^ Ad) 
valence scenario constraints (3.7) and (3.8) below, and which have to even vanish for the 
'unbroken' (Au = Ad) constraints (3.2) and (3.3). Therefore, for the time being, one has 
to resort to some, as far as possible, general and not too restrictive model assumptions 
concerning the breaking of the flavor-symmetry of the polarized light sea (8u ^ Sd), such 
as the proportionality (3.5). Only future polarized pp and pd Drell-Yan / u + pair and 
weak vector boson production experiments as well as polarized semi-inclusive DIS 
ep(d) — > e (tc, K)X experiments |5(| can provide us with direct measurements of the 
individual 5u(x, Q 2 ) and 5d(x, Q 2 ) distributions. 

We now have, instead of (3.2), 



Aq 3 = Au v (fi 2 ) - Ad v (fi 2 ) + 2[Aw(/i 2 ) - Ad(fi 2 )} =F + D 



(3.7) 



and on account of (3.3) and (3.6) 



Ag 8 (/i 2 ) 



Au v (fi 2 ) + Ad v (fi 2 ) + 2[Au(fi 2 ) + Ad(fi 2 )] 
3F-D + 2[Au(fi 2 ) + Ad(fi 2 )} 

as(m 2 ) 



(3.8) 
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where the imposed constraint (3.7) guarantees that the Bjorken sum rule (2.9) holds 
manifestly. Thus Eq. (1.5) becomes 



rr(Q 2 ) 



±^(F + D) + A(3F - D) + H (An(Q 2 ) + Ad(Q 2 )) 



1 - 



2\ 



7T 



(3.9) 

apart from a marginal contribution As(Q 2 ) = As(Q 2 ) < which is generated dynamically 
via the NLO evolution to Q 2 > /i^Lo even f° r the vanishing input in (3.6). Thus, in this 
case, only the total light-quark sea contribution in (3.9) has to be negative, Au(Q 2 ) + 
Ad(Q 2 ) < 0, in order to achieve the experimentally required reduction of the Ellis- Jaffe 
expectation (2.4). 

Our resulting input distributions can be parametrized as in (2.6) where now, instead of 
the unbroken Sq sea, we have a similar parametrization for 5u(x,[i 2 ) and 5d(x,[i 2 ) which 
are constrained by (3.5) ||51|| , together with the respective flavor-broken unpolarized input 
densities u(x, fj, 2 ) und d(x, n 2 ) taken from f2"5| . The parameters of our optimal LO densities 
at = 0.26 GeV 2 and the ones of the NLO(MS) densities at /4 LO = 0.40 GeV 2 are 
given in Table III. These optimal fits correspond to Xdf lo = 0.823 and Xdf nlo = 0.816 
similarly to the 'standard' scenario in Sec. 2. The quality of these LO and NLO fits 
to Af(x, Q 2 ) in the broken 'valence' scenario is practically identical to our new fit in 
the 'standard' scenario shown in Fig. 1 by the solid curves. The same holds true also for 
g\' n ' d {x, Q 2 ) shown in Figs. 2 und 3. The corresponding LO and NLO parton distributions 
at the respective input scales Q 2 = /i 2 G NLO in Eq. (2.6) fTj] with the 'valence' scenario 
fit parameters given in Table III are shown in Fig. 6, which are also compared with 



our reference unpolarized valence-like dynamical input densities of which satisfy the 
positivity constraint \5f\ < f. The polarized gluon densities turn out to be somewhat 
larger here, in particular in NLO, than the ones in the 'standard' scenario shown in Fig. 4. 
It should be furthermore emphasized that we always expect for the broken light-sea input 
densities to have a positive 5u and a negative Sd with \5d\ > 5u, i.e. 5u — 5d > and 
5u + 5d < 0. 
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In Fig. 7 we present the flavor asymmetry x(Su — Sd)(x, /i 2 ) separately, as obtained from 
Fig. 6, which compares favorably with predictions of the relativistic field theoretical chiral 
quark-soliton model |52| , |50| . Similar results have been obtained by a recent analysis ]53] 



based on the statistical parton model which are supposed to hold at a somewhat larger 
input scale Ql = M 2 ~ 0.9 GeV 2 . We note furthermore that the prediction 5g(x, M 2 ) = 
of this latter model is consistent with the results of the present analysis which do not 
exclude this possibility, cf. Fig. 3. The resulting NLO distributions at Q 2 = 5 GeV 2 are 
shown in Fig. 8 where they are also compared with our new NLO results obtained in the 
'standard' scenario as shown by the solid curves in Fig. 5. 

The first moments (total polarizations) Af(Q 2 ) of the polarized parton distributions 
of our fully flavor-broken 'valence' scenario and the resulting T^i^Q 2 ) are in LO given by 

Am = 0.664, Ad = -0.248, Au = 0.093, Ad = -0.261, As = As = 

Ag(fjl ) = 0.300, A#(5GeV 2 ) = 0.963, A#(10GeV 2 ) = 1.122 (3.10) 

which result in AS = 0.248 and 

r? = o.i4o, r? = -0.071. (3.11) 

Our NLO results are summarized in Table IV at some typical values of Q 2 . The nucleon's 
spin is carried almost entirely by the total helicities of quarks and gluons at the LO and 
NLO input scales 

iA£(/i 2 0;NL0 ) + As0u 2 o>NLO ) ~ 0.42, 0.48 (3.12) 

which is larger than the 'standard' scenario results (2.10), and thus a very small orbital 
contribution L q+g (fj,1 NLO ) ~ 0.08, 0.02 is required at the low input scales in order 
to comply with the sum rule (1.9). This is somewhat similar to our previous results 
[IJ, but again Ag(Q 2 ) is not strongly constrained by present data. Nevertheless it is 
intuitively appealing that this nonperturbative orbital (angular momentum) contribution 
to the helicity sum rule (1.9) vanishes at our low input scales, L q+g (fi 2 ) ~ 0. This is in 
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contrast to larger scales Q 2 > [i 2 where hard radiative effects give rise to sizeable orbital 
components due to the increasing kx of the partons, which eventually have to compensate 
in (1.9) the strongly increasing gluon polarization Ag(Q 2 ) ~ a~ 1 (Q 2 ) : in both scenarios 
we obtain, for example, Ag(10 GeV 2 ) ~ 1. 

Finally let us conclude with a few remarks concerning the flavor-symmetry breaking 
which was implemented in our broken 'valence' scenario via the entirely empirical relation 
(3.5). On rather general grounds one expects the product 

6q(x, /i 2 ) 6q(x, /i 2 ) = P(x) (3.13) 

to be a universal flavor-independent function P(x), since the effect of Pauli-blocking is 
only related to the spin (helicities) of quarks and antiquarks irrespective of their flavor 
degree of freedom. This implies Su(x, n 2 ) 5u(x, /i 2 ) = 5d(x, jj?) 5d(x, /i 2 ), i.e. Eq. (3.5). 
Furthermore, we have seen that the data select, within our 'valence' scenario with its 
totally flavor-broken polarized light sea densities in (3.5) and (3.6), the solution of Eq. 
(3.5) which satisfies P(x) > in (3.13) for q = u, d as can be seen in Fig. 6. This can be 
understood |50| as a consequence of the expected predominant pseudoscalar configuration 
29] , |54| of the quark-ant iquark pairs in the nucleon sea. In fact, the two relations uu ~ dd 



and Su 5u = 5d 5d at the input scale Q 2 = fi 2 can be rewritten as 

P p (x) = u + u + + U-U- ~ d + d + + d_d_ 

P a (x) = + u^u + ~ d + cL + d_d + (3-14) 

with P = P p — P a in (3.13) and the common helicity densities being given by Q ± = ( q 
± 8q )/2 where for brevity we have dropped the x-dependence. A predominant pseu- 
doscalar configuration of (gg)-pairs in the nucleon sea implies, via Pauli-blocking, that 
the aligned quark-quark configurations q + (q + q_) and g_(g_g + ) are suppressed relatively 
to the antialigned g + (g_g + ) and g_(g + g_) 'cloud' configurations, i.e. P p (x) > P a (x) which 
implies P(x) > in (3.13). The result for P p /P a , corresponding to our optimal fit, is 
shown in Fig. 9: clearly, this ratio will be maximal where xq(x,[i 2 ) and x5q(x,[i 2 ) are 
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maximal at x ~ 0.2 — 0.4, cf. Fig. 6, i.e. where the Pauli-blocking, Eq. (3.13), is most 
effective which is nicely exhibited in Fig. 9 in LO and NLO. 

It is interesting to mention that some of these expectations, which derive mainly 
from our light-sea flavor breaking relation (3.5), have been already confirmed by a recent 
entirely independent simultaneous analysis [TJ]] of polarized DIS and semi-inclusive deep 
inelastic scattering (SIDIS) asymmetry-data. In particular the more recent high precision 
SIDIS HERA-HERMES data [55 on h + production (h = rr,K dominantly) off a proton 
target, ep — > eh + X, seem to play a decisive role in favoring flavor-broken light sea 
densities 8u(x, Q 2 ) ^ d~d(x,Q 2 ), despite the fact that these asymmetry data on A\ v refer 
to rather small scales Q 2 ^ 1 GeV 2 . The reason for this discriminative power is, when 



combined with the data from inclusive DIS, due to the fact that A\ p is proportional Jl4 



besides to the dominant valence contribution, also to 5d — 45u, multiplied by a 'favored' 
fragmentation function, which is significantly more sensitive to 5u than to Sd. A clear 
preference for a positive Su has been observed ||14|| , which is very similar to our NLO 



Su shown in Fig. 6(b), and a flavor symmetric 'standard' light sea scenario seems to be 
strongly disfavored. 

We have calculated at NLO the spin asymmetries for semi-inclusive DIS using the well 
known theoretical SIDIS framework |56|, [14| together with our results for the polarized 
parton distributions of the 'standard' and 'valence' scenario with their flavor-symmetric 
and flavor-broken light sea densities, respectively, employing the fragmentation functions 
of |57| . (We did not use the alternative set of recent fragmentation functions suggested in 
|58|, since they refer to scales Q 2 larger than 2 GeV 2 .) The results for the relevant SIDIS 



asymmetry A\ are shown in Fig. 10. Although the high precision HERMES data [5^j 
seem to favor slightly the 'valence' scenario with its flavor-broken light sea, the results of 
the 'standard' scenario with its flavor-symmetric light sea cannot yet be ruled out. Both 
scenarios in Fig. 10 give rise to a comparable xV(9 data points) of 7.6 and 8.5 for the 
'valence' and 'standard' scenario, respectively. 
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4 Summary and conclusions 



All recent polarized DIS data, including the most recent SLAC-E155 proton data |21f . 
have been analyzed and studied within the 'standard' and 'valence' scenario in LO and 
NLO of QCD. The 'standard' scenario, characterized by (2.1) and (2.2), refers to the 
common simplified, but probably unrealistic, assumption of an SU(3)/ flavor-symmetric 
polarized light sea. The original 'valence' scenario characterized by (3.2) and (3.3), is 
now modified by employing a totally SU(3) f asymmetric polarized light sea 5u ^ 5d ^ 5s 
which leads to the modified constraints (3.7) and (3.8). Since inclusive polarized DIS data 
cannot fix the flavor-broken sea densities, we have modelled the flavor-asymmetric light 
sea densities 5u ^ 5d using a Pauli-blocking ansatz |3l| in (3.5), because similar 'blocking' 



effects can also explain the flavor-asymmetry of unpolarized sea densities (d > u). All 
our resulting polarized parton distributions respect the fundamental positivity constraints 
down to the low resolution scales Q 2 = fi1 Q = 0.26 GeV 2 and /Unlo = 0-40 GeV 2 . 
The polarized gluon distribution 5g(x, Q 2 ) is weakly constrained by present data in both 
scenarios. In particular, a vanishing gluon input 5g(x,fi 2 ) = is equally compatible 
with all present measurements of A±(x,Q 2 ) or g^(x,Q 2 ). Only a fully saturated (via 
the positivity constraint) gluon input 5g(x,fi 2 ) = ±g(x,fi 2 ) appears to be disfavored by 
present data. 

The presumably more realistic 'valence' scenario with its flavor-broken light sea quark 
distributions 5u ^ 5d (7^ 5s = 5s ~ 0) leads to a positive 5u(x, Q 2 ) density and a sizeably 
larger negative 5d(x, Q 2 ). These results are supported by a recent combined analysis [0 
of polarized DIS and semi-inclusive DIS data and agree with predictions of the relativistic 
field theoretical chiral quark-soliton model [50] and of the statistical parton model 



Present high statistics HERA-HERMES data [55 on semi-inclusive asymmetries 



for /i ± production off nucleon targets cannot, however, yet uniquely discriminate 
between our 'valence' scenario with flavor-broken polarized light sea densities and the 
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common 'standard' scenario with a flavor-symmetric light sea-quark distribution. 



A FORTRAN package containing our optimally fitted 'standard' and fully flavor- 
broken 'valence' NLO(MS) as well as LO distributions can be obtained by electronic mail. 



Acknowledgements 

We are grateful to Greg Mitchell for useful information concerning the recent E155 pro- 
ton data. This work has been supported in part by the 'Bundesministerium fur Bildung, 
Wissenschaft, Forschung und Technologie' and the 'Deutsche Forschungsgemeinschaft', 
Berlin/Bonn. W.V. is grateful to RIKEN, Brookhaven National Laboratory and the U.S. 
Department of Energy (contract number DE-AC02-98CH10886) for providing the facilities 
essential for the completion of his work. 



17 



References 

[1] M. Gliick, E. Reya, M. Stratmann, and W. Vogelsang, Phys. Rev. D53, 4775 (1996). 

[2] T. Gehrmann and W.J. Stirling, Phys. Rev. D53, 6100 (1996). 

[3] R.D. Ball, S. Forte, and G. Ridolfi, Phys. Lett. B378, 255 (1996); 

G. Altarelli, R.D. Ball, S. Forte, and G. Ridolfi, Nucl. Phys. B496, 337 (1997); Acta 
Phys. Pol. B29, 1145(1998). 

[4] K. Abe et al, SLAC-E154 Collab., Phys. Lett. B405, 180 (1997). 

[5] B. Adeva et al., CERN-SM Collab., Phys. Rev. D58, 112002 (1998). 

[6] C. Bourrely, F. Buccella, O. Pisanti, P. Santorelli, and J. Soffer, Prog. Theor. Phys. 
99, 1017 (1998). 

[7] D. de Florian, O.A. Sampayo, and R. Sassot, Phys. Rev. D57, 5803 (1998). 

[8] E.A. Leader, A.V. Sidorov, and D.B. Stamenov, Int. J. Mod. Phys. A13, 5573 (1998); 
Phys. Rev. D58, 114028 (1998); Phys. Lett. B445, 232 (1998); B462, 189 (1999); 
Phys. Lett. B488, 283 (2000). 

[9] L.E. Gordon, M. Goshtasbpour, and G.P Ramsey, Phys. Rev. D58, 094017 (1998). 

[10] S. Incerti and F. Sabatie, SLAC-E154/155 Collab., CEA-Saclay PCCF-RI-9901 
(1999). 

[11] S. Tatur, J. Bartelski, and M. Kurzela, Acta. Phys. Pol. B31, 647 (2000); 
J. Bartelski and S. Tatur, |hep-ph/00042511 (2000). 



[12] D.K. Ghosh, S. Gupta, and D. Indumathi, Phys. Rev. D62, 094012 (2000). 
[13] Y. Goto et al., Asymmetry Analysis Collab., Phys. Rev. D62, 034017 (2000). 
[14] D. de Florian and R. Sassot, Phys. Rev. D62, 094025 (2000). 

18 



[15] SLAC-Yale E80 Collaboration, M.J. Alguard et al., Phys. Rev. Lett. 37, 1261 (1976); 
G. Baum et al., ibid. 45, 2000 (1980); 

SLAC-Yale-E130 Collaboration, G. Baum et al, Phys. Rev. Lett. 51, 1135 (1983). 

[16] CERN-EM Collaboration, J. Ashman et al., Phys. Lett. B206, 364 (1988); Nucl. 
Phys. B328, 1 (1989). 

[17] CERN-SM Collaboration, B. Adeva et al, Phys. Rev. D58, 112001 (1998). 

[18] SLAC-E143 Collaboration, K. Abe et al., Phys. Rev. D58, 120003 (1998). 

[19] SLAC-E142 Collaboration, P.L. Anthony et al., Phys. Rev. D54, 6620 (1996). 

[20] DESY-HERMES Collaboration, K. Ackerstaff et al., Phys. Lett. B404, 383 (1997). 

[21] DESY-HERMES Collaboration, A. Airapetian et al, Phys. Lett. B442, 484 (1998). 

[22] SLAC-E154 Collaboration, K. Abe et al., Phys. Rev. Lett. 79, 26 (1997). 

[23] SLAC-E155 Collaboration, P.L. Anthony et al., Phys. Lett. B463, 339 (1999). 

[24] SLAC-E155 Collaboration, P.L. Anthony et al, SLAC-PUB-7994, July 2000 T/E 
( |hep-ph/0007248| ). 



[25] M. Gluck, E. Reya, and A. Vogt, Eur. Phys. J. C5, 461 (1998). 

[26] A.D. Martin, R.G. Roberts, W.J. Stirling, and R.S. Thorne, Eur. Phys. J. C4, 463 
(1998); Eur. Phys. J. C14, 133 (2000). 

[27] H.L. Lai et al., CTEQ5 Collab., Eur. Phys. J. C12, 375 (2000). 

[28] H. Botje, Eur. Phys. J. C14, 285 (2000). 

[29] For recent reviews see, for example, S. Kumano, Phys. Rep. 303, 183 (1998); 
J. Speth and A.W. Thomas, Adv. Nucl. Phys. 24, 83 (1998). 



19 



[30] M. Gliick and E. Reya, Mod. Phys. Lett. A15, 883 (2000) ( |hep-ph/0002T82|) . 

[31] W. Furmanski and R. Petronzio, Z. Phys. Cll, 293 (1982). 

[32] M. Gliick, E. Reya, and W. Vogelsang, Nucl. Phys. B351, 579 (1991). 

[33] Particle Data Group, C. Caso et al, Eur. Phys. J. C3, 1 (1998). 

[34] M. Karliner and H.J. Lipkin, Phys. Lett. B461, 280 (1999). 

[35] J. Ellis and R.L. Jaffe, Phys. Rev. D9, 1444 (1974); 10, 1669 (E) (1974); 
M. Gourdin, Nucl. Phys. B38, 418 (1972). 

[36] S.J. Brodsky and I. Schmidt, Phys. Lett. B234, 144 (1990); S.J. Brodsky, 
M. Burkardt, and I. Schmidt, Nucl. Phys. B441, 197 (1995); 
F.E. Close and D. Sivers, Phys. Rev. Lett. 39, 1116 (1977). 

[37] Note that our ansatz violates (slightly) the 'helicity retention' concept developed 
in Ref. ||36[| , by which Sf(x,Q 2 )/f(x,Q 2 ) — > 1 for x — > 1. Such relations follow in 
perturbative QCD and need not necessarily be fulfilled by our non-perturbative input 
at Q 2 = fi 2 . In any case, we could always appropriately modify our ansatz at very 
high x so as to satisfy the 'helicity retention' conditions. Enforcing in our fits all f3f 
to vanish in Eq. (2.6) increases the total x 2 by several units, giving rise to worse fits 
than our optimal ones. 

[38] We have refrained from attempting to extract information from measurements of 
the spin asymmetry for photoproduction of oppositely charged hadron-pairs pre- 
sented by the HERMES collaboration [A. Airapetian et al., Phys. Rev. Lett. 84, 
2584 (2000)]. Such data could in principle constrain the spin gluon density Sg(x, Q 2 )] 
however, the smallness of the hadron transverse momenta accessed in this experi- 
ment makes an interpretation of the data in terms of QCD hard scattering difficult. 
In fact, as stated by the HERMES collaboration, their published gluon asymmetry 



20 



6g(x, Q 2 )/g(x, Q 2 ) = 0.41 ± 0.18 at x ~ 0.17 and Q 2 = (p 2 T ) = 2.1 GeV 2 , which im- 
plies x5g(x, Q 2 ) ~ 0.35±0.15, actually compatible with our results in Fig. 5, does not 
include systematic theoretical errors. These could be sizeable in the low transverse 
momentum region considered. 

[39] J.D. Bjorken, Phys. Rev. 148, 1467 (1966); ibid. Dl, 1376 (1970). 

[40] H.Y. Cheng, Int. J. Mod. Phys. All (1996) 5109; Phys. Lett. B427 (1998) 371. 
The explicit transformation of the splitting functions (anomalous dimensions) and 
coefficient functions from the MS to the CI (JET) scheme can be found in this latter 
reference as well as in the third reference of [8]. 

[41] D. Miiller and O.V. Teryaev, Phys. Rev. D56, 2607 (1997). 

[42] H. Abramowicz et al., CDHSW Collab., Z. Phys. C15, 19 (1982). 

[43] CCFR Collab., C. Foudas et al., Phys. Rev. Lett. 64, 1207 (1990); 

SA. Rabinowitz et al., ibid. 70, 134 (1993); A.O. Bazarko et al., Z. Phys. C65, 189 
(1995). 

[44] H.J. Lipkin, Phys. Lett. B256, 284 (1991); 337, 157 (1994); 
J. Lichtenstadt and H.J. Lipkin, ibid. 353, 119 (1995). 

[45] A. Baldit et al., NA51 Collab., Phys. Lett. B332, 244 (1994); 

EA. Hawker et al, E866 Collab., Phys. Rev. Lett. 80, 3715 (1998). 

[46] CA. Miller, HERMES Collab., Proceedings of DIS '99, Nucl Phys. B (Proc. Suppl.) 
79, 146 (1999). 

[47] R.D. Field and R.P Feynman, Phys. Rev. D15, 2590 (1977). 

[48] S. Kumano and M. Miyama, Phys. Lett. B479, 149 (2000); 

B. Dressier et al., Univ. Bochum RUB-TPII-13/99 flhep-ph/ 9910464ft ; 



21 



C. Bourrely and J. Soffer, Phys. Lett. B314, 132 (1993); Nucl. Phys. B423, 329 
(1994). 

[49] T. Morii and T. Yamanishi, Phys. Rev. D61, 057501 (2000). 
[50] B. Dressier et al, Eur. Phys. J. C14, 147 (2000). 

[51] The easiest and consistent way to implement directly (even for Mellin-n moments) 
the constraint (3.5) in a fit routine is to fit, for example, 8u, 5d and Su according to 
the ansatz (2.6) which fix 5d(x,fi 2 ) via (3.5) and (3.4): 5d/d ~ (Su/u) (5u/u)/(8d/d) 
where the unpolarized densities are taken from Ref. P5| . Thus N d = N U N U /N d , 
aj = a u + a a - a d and #j = f3 u + f3 u - (3 d . 

[52] D. Diakonov et al., Nucl. Phys. B480, 341 (1996); Phys. Rev. D56, 4069 (1997); 
M. Wakamatsu and T. Kubota, Phys. Rev. D60, 034020 (1999); 
K. Goeke et al, Univ. Bochum RUB-TPII-18/99 ( |hep-ph/000T272l) . 

[53] R.S. Bhalerao, Phys. Lett. B380, 1 (1996); 387, 881 (E) (1996); talk given at the 
Int. Conf. on QNP 2000, Adelaide, Feb. 2000, Nucl. Phys. A680, 62 (2000); 
R.S. Bhalerao, N.G. Kelkar, and B. Ram, Phys. Lett. B476, 285 (2000). 

[54] N.N. Nikolaev et al., Phys. Rev. D60, 014004 (1999). 

[55] K. Ackerstaff et al., HERMES Collab., Phys. Lett. B464, 123 (1999). 

[56] D. de Florian, C.A. Garcia Canal and R. Sassot, Nucl. Phys. B470, 195 (1996). 

[57] S. Kretzer, Phys. Rev. D62, 054001 (2000). 

[58] B.A. Kniehl, G. Kramer and B. Potter, Nucl. Phys. B582, 514 (2000). 



22 



Table I. The parameters of the LO and NLO input parton distributions, as defined 
in Eq. (2.6), at filo = 0-26 GeV 2 and /i^LO = 0-40 GeV 2 , respectively, for the 'standard' 
scenario as obtained from fits to the data in [15 - 24]. 





standard scenario 




LO 


NLO 


AT 


n 8^1 


1 m q 


U. u 


0.45 


0.52 


Pu 





0.12 


N d 


-0.734 


-0.669 


ad 


0.49 


0.43 


(3 d 


0.03 







-0.587 


-0.272 




0.68 


0.38 


fa 








N s 


1 


1 


N g 


1.669 


1.419 


a g 


1.79 


1.43 


P 9 


0.15 


0.15 


X 2 /209 data pts. 


174.9 


169.8 



Table II. First moments (total polarizations) Af of polarized NLO parton densities 
Sf(x, Q 2 ) and gi' n (x, Q 2 ), defined in (1.3) and (1.4), as obtained in the 'standard' scenario. 
The marginal differences between Au and Ad at Q 2 > fi 2 , generated dynamically by the 
NLO evolution, are not displayed. 



Q 2 (GeV 2 ) 


Au 


Ad 


Ag 


Ag 
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rn 
1 1 


f^NLO 
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-0.054 


1 


0.861 
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0.132 


-0.062 


10 


0.859 


-0.406 


-0.064 


0.828 


0.197 


0.133 


-0.063 
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Table III. The parameters of the LO and NLO input distributions, as defined in 
Eq. (2.6), with q to be identified with u, at nl = 0.26 GeV 2 and /4 LO = 0.40 GeV 2 , 
respectively, for the 'valence' scenario. We have fitted the broken light-sea input density 
5u(x,fi 2 ) and fixed 5d(x,[i 2 ) via (3.5) and (3.4) |51|j . 
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Table IV. First moments (total polarizations) Af of polarized parton densities 
Sf(x, Q 2 ) and gi' n (x, Q 2 ), defined in (1.3) and (1.4), as obtained in the fully flavor-broken 
'valence' scenario. The marginal finite As = As(Q 2 > fJ^ LO ) are generated dynamically 
by the NLO evolution. 
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Figure Captions 



Fig. 1. Comparison of our NLO results for A^(x,Q 2 ) as obtained from the fitted input at 
Q 2 = /^nlo f° r the 'standard' SU(3)/ flavor-symmetric sea scenario [Eq. (2.6) and 
Table I] with present data [16 - 24] . The Q 2 values adopted here correspond to the 
different values quoted in [16 - 24] for each data point starting at Q 2 > 1 GeV 2 at 
the lowest available x bin. Our old NLO GRSV95 fit [[1J is shown for comparison as 
well (dashed curves). Our present LO fit is very similar to the NLO one shown by 
the solid curves. 

Fig. 2. Comparing the Q 2 -dependence of our LO and NLO 'standard' scenario results with 
recent data [16 - 18, 21, 24] on g[(x,Q 2 ). To ease the graphical representation we 
have multiplied all results at the various fixed values of x by the numbers indicated 
in parentheses. 

Fig. 3. The x-dependence of at Q 2 = 5 GeV 2 in the NLO 'standard' scenario. Different 
choices of the gluon input 5g(x, /^nixO m (2-6) are shown by the dashed and dotted 
curves as indicated. Allowing our optimal total x 2 to change by one unit, 5x 2 = ±1, 
results in the shaded areas shown. The data are taken from [17, 18, 20, 22 - 24]. 

Fig. 4. (a) Comparison of our fitted 'standard' LO input distributions in Eq. (2.6) and 
Table I at /x£ Q = 0.26 GeV 2 with our previous old GRSV95 fit g] and with the 



unpolarized dynamical GRV98 input densities of |25]. (b) The same as in (a) but 
for the NLO input densities at ^ LO = 0.40 GeV 2 . 

Fig. 5. The polarized LO and NLO(MS) distributions at Q 2 = 5 GeV 2 in the 'standard' 
scenario, as obtained from the input densities at Q 2 = /i 2 NLO in Fig. 4. Our old 
NLO results [0 are shown for comparison. 

Fig. 6. (a) Comparison of our fitted fully broken 'valence' scenario LO input distributions 
in Eq. (2.6) and Table III at /z£ = 0.26 GeV 2 with the unpolarized GRV98 
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input densities of (b) The same as in (a) but for the NLO input densities at 
/4 LO = 0.40 GeV 2 . 

Fig. 7. LO and NLO results for the difference of the broken light-sea input densities 
5u(x,fi 2 ) and 5d(x,/i 2 ) in the 'valence' scenario as obtained from Fig. 6. The pre- 



diction of the chiral quark-soliton model is taken from K. Goeke et al. [52]. 



Fig. 8. The polarized LO and NLO(MS) distributions at Q 2 = 5 GeV 2 in the fully flavor- 
broken 'valence' scenario, as obtained from the input densities at Q 2 = /^lo, nlo 
in Fig. 6. For comparison the new NLO results of the 'standard' (unbroken sea) 
scenario are shown as well by the dashed curves (which coincide with the solid curves 
of Fig. 5.). 

Fig. 9. The ratio of P p {x) and P a (x), defined in (3.14), at the LO and NLO input scales 
yU 2 NLO which demonstrates the Pauli-blocking of the disfavored antiparallel q±q T 
configurations relative to the favored parallel q±q± configurations as discussed in 
the text. 

Fig. 10. NLO predictions for the semi-inclusive DIS asymmetry A^l for h + production off a 
proton target within the 'valence' and 'standard' scenario. The HERMES data are 



taken from |55|. 
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